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We study effective actions for order parameter fluctuations at low temperature in 
layered d-wave superconductors such as the cuprates. The order parameter lives on 
the bonds of a square lattice and has two amplitude and two phase modes associated 
with it. The low frequency spectral weights for amplitude and relative phase fluc- 
tuations is determined and found to be subdominant to quasiparticle contributions. 
The Goldstone phase mode and its coupling to density fluctuations in charged sys- 
tems is treated in a gauge- invariant manner. The Gaussian phase action is used to 
study both the c-axis Josephson plasmon and the more conventional in-plane plas- 
mon in the cuprates. We go beyond the Gaussian theory by deriving a coarse-grained 
quantum XY model, which incorporates important cutoff effects overlooked in previ- 
ous studies. A variational analysis of this effective model shows that in the cuprates, 
quantum effects of phase fluctuations are important in reducing the zero temperature 
superfluid stiffness, but thermal effects are small for T<T C . 

PACS numbers: 74.20.-z, 74.20.De, 74.72.-h, 74.25.Nf 



I. INTRODUCTION 



The high temperature cuprate superconductors (SCs) differ from conventional SCs in several respects: a d-wave gap 
with gapless quasiparticle excitations, a small superfluid phase stiffness, a short coherence length and strong electron 
interactions. It is therefore of interest to examine some of these unconventional aspects and their interplay in simple 
models of high T c systems. With this motivation, we study in this paper, the low temperature collective properties 
of charged, layered d-wave SCs with a short coherence length and small superfluid stiffness. 

The superfluid phase stiffness, D s — n s /m*, is a fundamental property characterizing SCs which is directly 
related to A, the experimentally measured magnetic penetration depth Q in the London limit . The low temperature 
behavior of D s contains information on the low-lying excitations in these systems. Experimentally, A(T) is found to 
increase linearly with T in the high T c SCs implying a linearly decreasing D s . This linear drop in D s has been 
attributed to quasiparticle excitations near the nodes of the ci-wave gap. Alternatively, it has been suggested that 
this effect could arise entirely from classical thermal phase fluctuations ||,U and quasiparticles can be ignored 
It is then clearly of interest to identify the important low energy excitations in these systems, from the point of view 
of understanding the penetration depth data, as well as other thermodynamic properties and response functions. 

From a theoretical perspective, the physics of a system with a small superfluid stiffness and short coherence length 
has been studied in detail in case of neutral s-wave SCs HQ. In this case, the fermionic excitations and fluctuations 
in the order parameter amplitude are gapped, and phase fluctuations are the only important excitation at low tem- 
perature. It is of interest to compare this with the behavior in models which support an anisotropic order parameter 
with low lying fermionic excitations, such as a ci-wave SC. 

We approach the problem by deriving and analyzing effective actions for a d-wave SC within a functional integral 
framework, which allows us to focus on the collective (order parameter) degrees of freedom. We note that our 
effective actions are derived by looking at fluctuations around a BCS mean field solution. We believe that such an 
approach is valid for the SC state of the high Tc materials, at least for T <C T c . There is considerable experimental 
evidence for sharply defined quasiparticle excitations about the d-wave SC state, and thus the ground state and low- 
lying excitations appear to be adiabatically connected to their BCS counterparts. We thus expect that while strong 
correlations will modify the coefficients of the phase action, they will not change its qualitative form. 

Our main results can be summarized as follows: 

1. We find that the d-wave SC state is characterized in terms of an order parameter which lives on the bonds 
of a square lattice. The bond order parameter leads to two amplitude and two phase modes in contrast to s-wave 
SCs. One of the phase modes is identified as the usual (Goldstone) phase mode while the other, which we call the 
"bond phase" , is the relative phase between the x and y bonds at a site. The latter can be thought of as representing 
fluctuations from the d-wave state towards an extended s-wave state. The amplitude and bond-phase fields have spin 
zero and couple to the particle-particle channel. 
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2. We study the spectral weight for fluctuations of the amplitude and bond-phase fields and find that they are 
not gapped but rather exhibit power laws down to zero energy. However, the low energy spectral weight in these 
fluctuations is very small compared to the quasiparticle contribution. 

3. We derive an effective Gaussian action for the usual phase variable in charged systems, since this couples 
directly to the electromagnetic potentials. The large in-plane plasma frequency, which is relatively unaffected by 
superconductivity, and the low energy c— axis Josephson plasmon at zero and finite temperatures are studied in a 
unified manner within the same formalism. We emphasize the relation between unusual aspects of the c-axis optical 
conductivity and the Josephson plasmon. We also discuss the plasmon dispersion in layered systems. 

4. We extend the above formalism to consider the effect of phase fluctuations beyond Gaussian level on the superfiuid 
stiffness in charged systems. The quantum XY model || is usually used for such an analysis, motivated by studies 
of Josephson junction arrays and granular SCs. We emphasize that there are important differences when considering 
low temperature bulk SCs, and derive a quantum XY phase action suitable for our problem, correctly taking into 
account appropriate momentum and frequency cutoffs, missed in earlier studies. 

5. The low temperature renormalization of D s by phase fluctuations is studied within a self-consistent harmonic 
approximation. For parameter values relevant to the cuprates near optimal doping, quantum phase fluctuations are 
shown to lead to a sizeable renormalization of the superfiuid stiffness. However, thermal fluctuations are found to 
have no effect at low temperatures, unlike the results of earlier studies These studies focussed on the effect of 
thermal phase fluctuations, but Coulomb effects were considered to be unimportant, in contrast to the present work. 

6. As part of our analysis, we also touch upon certain formal issues which may be of some general interest. Among 
these are: (a) how gauge invariance can be understood in a simple manner within the functional integral language; 
(b) the role of the linear time derivative term ipdr9 in the phase action; and (c) the problems involved in deriving 
local phase actions which respect 2ir periodicity starting from a fermionic model. 

The paper is organized as follows. In Section II, we present the Hamiltonian for our model, and discuss the effective 
action and mean field theory in Section III. Section IV contains a discussion of fluctuations of the amplitude and 
bond-phase fields with some of the details discussed in Appendix A. In Section V, we turn to phase fluctuations and 
derive effective phase-only actions for neutral and charged systems. The linear time derivative term in the action 
which arises in this context is briefly discussed in Appendix B. We then derive gauge invariant density and current 
correlations, leaving details of the algebra to Appendix C. In Section VI we discuss collective in-plane and c— axis 
plasmons. In Section VII we present the derivation of a quantum XY model appropriate for charged, layered SCs. 
We analyze this action and compute the renormalization of the phase stiffness by longitudinal phase fluctuations in 
Section VIII and discuss experimental implications. We conclude in Section IX with a discussion and summary of our 
results. 



II. THE HAMILTONIAN 

We consider a system of fermions with kinetic energy K = J2k a & c k <j c K<? (where £k = £k — M with the 2D 
dispersion and /i the chemical potential) interacting via a separable potential which is attractive in the d-wave channel. 
We will show that in coordinate space this interaction leads to the superexchange term of the t-J model. 

Let us begin with 

tfpair = -jfJ2 ^( k )^rf( k ')c k + q /2T C -k+q/2i C -k'+q/21 C k'+q/2T (!) 
k,k',q 

where <Pd(k) = (cosfe^ — cosfcj,). and we work on a 2D square lattice with lattice spacing a. Q denotes the volume 
of the system with N sites. We set a = 1 in most equations, but retain it in some for the sake of clarity. Fourier 
transforming to real space, we get 

H pail = -| B t,r' B r,r' = | (S r • S r > - -n r rv). (2) 

(r,r') <r,r') 

(The prime on _ff pa i r is omitted in going from ([!]) to (J2J) for reasons explained below). Here (r, r') are nearest neighbour 
sites, and B r r , = c\ ^c T , ^ — c r ^c\, ^ creates a singlet on the bond (r, r'), while S r and n r are the spin and number 
operators. Clearly this is the interaction term of the t — J model with g = 2 J. 

There is a subtlety involved here; on transforming (g) back to k— space we do not recover the original expression 
(|l|) we had started with. Instead we obtain, using from now on g = 2 J, 
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#pair = --^ [^( k )^( k ') + ^( k )¥ , s( k ')]4+q/2T C -k+q/21 C -k'+q/2lCk'+q/2T (3) 

k,k' : q 

where <^ s (k) = (cosfc^ + cosk y ). The reason for two different k-space interactions leading to the same real space 
expression is the following operator identity on a 2D square lattice: 

E [tp d (k)tp d (k') - if s (k)ip s (k')} cj c +q/2T C -k+q/2i C - k '+q/ 2 i Ck '+ < l/ 2 T = ( 4 ) 
k,k',q 

It can be shown that (g) and <j3j> both lead to the same self consistent BCS gap equation. Thus we will use the 
interaction in (^), and not ([!]). From the form of (^|), it is clear that H pa i r has attraction in both the d-wave channel, 
with a fd(k) order parameter, and in the extended s-wave (s*) channel, with a <p s (k) order parameter. 

We will now analyze the Hamiltonian H = K + impair- Later, we will also add to it the Coulomb interaction 
appropriate to layered systems (see Section V.B). 



III. MEAN FIELD THEORY 



The partition function at a temperature T is written as the standard coherent state path integral with the ac- 
' ' ' ^ 1 1 " ! i ! We decouple H pail with a complex field A r y(r) using the Hubbard- 



tion J Q dr Y^,* c l,ad T c r ,a + H (c, c' 
Stratonovich transformation: 



exp(^< r ,(r)B r , r ,(T)) = J 23(A A*) exp [-£(r, r'; r)] 



(5) 



where 



We thus obtain the action 



£(r,r';r) = ^|A r .^(r)| 2 - \ ( A T , r , (tJB* (r) + h.c 



8J 



(6) 



S = 



l/T 



(It 



E <i» & + a) c k , CT (r) + E r '; ^) 



k.rr 



<r,r'> 



(7) 



The fermion fields can then be integrated out to obtain the effective action exp (— S c g [A, A*]) = J V (c, c^) e . 

The (i-wave saddle point is given by A r r+ x(r) = — A rjr _|_y(r) = A^, a (r, r)-independent real number, obtained 
from SS c ff/6Ad — 0, which leads to the BCS gap equation 



J N ^ 



yg( k ) 

2E k 



tanh(S k /2T) 



(8) 



where — yS~+ A 2 , and Ak = Ad(/3d(k)/2. The same gap equation can be obtained by starting from the momentum 
space potential in (H) and decoupling in the (i-wave channel. 



Given H r , 




one could equally well look for possible extended s-wave (s*) saddle points with A r r+ x(r) 



A r r+ y(r) = A s . However, for our choice of dispersion |9| (which includes nearest and next nearest neighbour hopping 
with opposite signs) we have found by numerical solution of the gap equation that the d-wave saddle point is stable 
relative to the s* solution. The reason for this can be seen as follows: v? s (k) is small over most of the Fermi surface 
for the fillings of interest, while fd(k) vanishes only on the nodes. Thus the condensation energy gained by the s* 
state is smaller than the d-wave state. Further, if we consider the large on-site repulsion between electrons (which we 
have not done here, but is certainly an essential part of the full t — J model) and demand (c^cj , ) = 0, then we have 
^ k Ak/2i?k = 0. This is automatically satisfied for the d-wave state at any filling, but not for the s* state. In this 
work, we rely on the former "Fermi surface effect" to stabilize the d-wave state. 
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IV. FLUCTUATIONS 



To treat fluctuations in the order parameter we write A rir <(r) = |A r . I .'(r)|e t * r . r '( T ). The phase $ r ,r+x( T ) = and 
^r.r+y (t) = 7t at the ci-wave saddle point. We now divide the phase field into two parts; following Ref |1C[| , we set 
$r,r+£'(V) = 0(r, r) and $ r , r+ y(r) = tt + 0(r, r) + 0(r, t). 

We next assume that the spatial variation of 9(r, t) is small on the scale of the lattice spacing, which allows us 
to set $ r , r+ i(r) « i [0(r,r) + 0(r + x,r)] and $ P)P+fi (r) s» 7r + </>(r, r) + ± [0(r,r) + 0(r + y, r)]. While we lose the 
0(r, r) — > 0(r, t) + 27r invariance of the action with this approximation, it is nevertheless useful in isolating that part 
of the phase field which couples to electromagnetic fields as we will see below. We can now transform to new fermion 
variables given by cJ(t) — ► c r (r)e~' ie ( r ' r )/ 2 [[IT 12|. As a result of this "gauge transformation" the action of (jl]) gets 
modified to 



S : 



,1/T 

/ dr [Co 
Jo 



A] 



(9) 



with 



-i8(r,r)/2 



(d T - fi)c r Ary 0{r - T)/2 r ') \4Ar)c rl Ar)e-^- T) - e ^'- T)]/2 + h.c.] (10) 



where t(r — r') is the hopping matrix element between point r and r', so that = ^2 r _ r , e lk ' r r H(r — r') and 
£l = ^7 E l A ^' Ml 2 - J E l A r,r+x(r)| « p+ *(t) + h.c.) + \ ]T |A p , p+ ^t)| (i?t ^ (r) e W) + h.c.) . (11) 

(r.r') r r 

In the following Sections we shall integrate out the fermions and examine the resulting effective actions for the 
amplitude, (f> and 9 fields. 



A. Amplitude fluctuations 

Amplitude fluctuations can be considered by setting |A PiP _|_&(t)| = A<j(l + i] a (r,T)) in (|il"|), where a — x,y. The 
transformation from A Q , A* to rj a ,9,(f> has a Jacobian 4A^(l+r/ x )(l+ryj,) at every point (r, r), leading to an additional 
term 



1 jT 1 jT 1 1 T 

[ dr C 2 =T [ dr Vln(l + ?7 a (r,T)) wT / dr V 
Jo Jo _ „ Jo „, 



(12) 



in the action in (||). For (q, cj) ^ (0,0), the linear term in ( |l2| ) can be set to zero and only the quadratic term 
contributes. However, even this term is zero at T — and can be ignored at low T. 

From (|lTJ) we see that the spin zero amplitude fields r\ a couple to singlet pairs. Their coupling to the </>- field can 
be shown to be small at small momentum and frequency. In particular, for static uniform distortions of <j> and r] a the 
energy has to be even under cj> — > — and terms like 4>rj a cannot appear in the action on integrating out the fermions. 
The mixing of rj a with the phase 9 and electromagnetic potentials can also be shown to be negligible at small q, uj 
since rj a couples to the particle-particle channel while the 9 and electromagnetic potentials couple to the particle-hole 
channel. The mixing then involves integrals over products of ordinary and anomalous Green functions which vanish 
using particle-hole symmetry near the Fermi surface and the k— dependence of At- This lack of mixing of amplitude 
and phase modes is similar to the well known weak coupling result for s-wave superconductors J7| . 

Unlike the s-wave case, however, the amplitude excitations have low frequency spectral weight in c£-wave SCs which 
we now estimate. Setting the phase fields to saddle point values and integrating out the fermions we obtain an 
effective action for amplitude fluctuations. Transforming to new variables -q s = (rj x + rj v )/y/2 and r\d = (Vx ~ Vy)/V^, 
diagonalizes the action for rj fields at q = and is a good starting point to consider small q fluctuations. We obtain, 
to Gaussian order, S[rf\ — ^ 5Z q „ i— s d Vi-M^ (<l, iuj n )r]i, where ui n — 2-kuT. We have made the approximation of 
ignoring the coupling between 7/ s and which can be shown to be negligible at small momentum. 

The low energy density of states is given by Ni(u) — jj ^ q ImMi(q, w)/-k with |q a |, \<\ y \ < 7r/£ - The cutoffs arise 
since the fluctuations must have an energy lower than the condensation energy as discussed in more detail later, in the 
context of phase fluctuations. From a numerical calculation of N St d(co) we find that N s (oj)/N qp (u>) ~ uj 2 /(A d v F /a) 



4 



and Nd(cu)/N qp (iv) ~ uj 4 /(A^v F /a) where N qp (uj) = k F tu / {irv F Ad) is the density of states per spin for quasiparticle 
excitations. These results can also be understood from an approximate analysis of the form of M S)( 2(q, u) discussed 
in Appendix A. 

We see that both N s (ui) and Nd(u>) are much smaller than N qp (u>) for uj -C Ad,and thus conclude that amplitude 
fluctuations are unimportant for low temperature properties which will be dominated by the quasiparticle contribution. 



B. The "bond-phase" field <f> 



We next study the field <fi. We note from (|Tl|) that a uniform <fi = n would lead to extended s-wave (s*) order. One 
can therefore think of (f> as representing fluctuations of s* character about the d-wave saddle point. From ( pT| ) we see 
4> is a spin zero field which couples to pairs. 

For reasons similar to those explained above for amplitude fluctuations, the coupling of tf> to other fields is 
weak and may be ignored at low momentum and frequency. We therefore derive a Gaussian action for (f> by 
setting 9 = rj s =0, their saddle point values, and integrating out the fermions. This leads to the action 
S[(f>] = y ^2 n 1 (q, w„)|0(q, uj n )\ 2 . Since the field has low frequency spectral weight, we compute its density 
of states N^cu) = Yl\q\<w/£ -^ rn -^>( t l> w + i0 + ). From numerical calculations, as well as simpler approximations 
discussed in Appendix A, we find N ^(uj) / N qp {uj) ~ u> 2 / '(A4V F / 'a) . We thus see that <f> fluctuations are much less 
important than quasiparticles at low temperatures. 



V. PHASE FLUCTUATIONS 



From action of (|10| ) we see that uniform shifts in 9 do not cost any energy, and 9 is the Goldstone mode of the 
superconducting state. 

We now obtain the action for 9~ fluctuations coupled to fermions, setting r) a = <j) = (their saddle point values) in 
( |l0| ) and (pi]). For slow spatial fluctuations, the deviation from the mean field action, obtained from (|l0|) with 9 = 0, 
is given by 

6S[J, c, 0] = ± [ 1/T dr 4,»Ck-q,*M [idr - i(& - &-,)] q (r) (13) 

1/T 

dT Y i lff Wck-q-q', ff (T)^(T)^(T) & + 6,-q-q' - &-q - ^k-q') 
k,q,q',CT 

Using 9q(r + 1/T) = 0q(r), an assumption discussed in detail in Appendix B, and making a small q expansion, we 
arrive at 

ss[c\ c , 9} = c L c *-«.A K - iv kQ q Q ] - Y c L c fc-«-9',^ m ^( k )q«q/3 6l ? 6, 9' ( 14 ) 

k,q,<7 k.q.q' .{J 

where k = (k, iv m ) and q = (q, iu> n ) with v m = (2m + 1)ttT and u> n = ItvkT . We have used £k-q — £k — Vk a qa + 
k m ap0 i ) ( la < ll3 + • ■ ■ where Vk Q = d^k/dk a is the velocity and m~^(k) = d 2 ^/dk a dkp is the inverse mass tensor. 




A. Neutral Systems 



For neutral systems we integrate out the fermions in (p_4[), using a cumulant expansion |13| controlled by small 
spatial and temporal gradients in 9, leading to 

Sneutral [9] = ^ £ £ [- X 0^ + Ao^qV] <?(q, t« n )0(-q, -<W n ). (15) 

Here, xo = — (pIq? * w n)p( — Qj — i^ n )) is the mean field density-density correlation function given by 



5 



Xo(q, iu> n ) = ^ - / - /')(««' + ««') 



zw„ - E - E' iu> n + E + E' 



^^(/-/OK-m') 



iw„ - £ + iw„ + E - E' 



(16) 



and A< 



a/3 _ 1 



nEk^WW ~ ^(j a (qj^n)j (~qi"* Li; n))o is the mean field phase stiffness, with (ra k ) 



(1 — ^k/£ ; ktanh(_Ek/2T)) and the paramagnetic current correlator 

™ (M3>0 = 75 XI V k« V k/3(l - / - /' - W') 



^v kQ v k(3 (/' - f)(vv' + uu') 



iuj n — E — E' iu n + E + E' 
vv' uu 1 



E-E> iu„ -E + E' 



(17) 



E, u and v refer to standard BCS notation, / = f(E) is the Fermi function, E 1 = i? k - q with E = E^ and similarly 
for other primed variables. Analytically continuing iu> n — ► u> + i?7, and working at T = |l4[ , we obtain in the limit 
q — > 0, a; — > 0, the mean field superfluid stiffness 



A o Q/3 (T = 0)- 



'a/3 



1 



in 



D s °(T = 0)<^ 



and the mean field compressibility 



*o(r = o) = -i£~ k 



B2 



(18) 



(19) 



We note that the effective action (JlJ) is appropriate for phase distortions whose energy is smaller than the conden- 
sation energy -E con d = \D q s {-k / '£ ) 2 ■ If this energy (density) is exceeded the d-wave BCS saddle point would become 



= E, 



cond ■ 



In the 



unstable. This leads to the following restrictions in (|15|): |q| < q c = 7r/£o an d <-<J n < w c with ku 
BCS limit, i? con d — A^/v F k F which translates into v F q c ~ uj c ~ A^. 

We emphasize that the coefficients in the phase action are not the physical correlation functions. The gauge invariant 
correlation functions, obtained by including the effect of Gaussian phase fluctuations, are given by 



and 



A a/3 = A aft 



Ao^Ao^q^q" 



q Q q /3 A 



a/3 



Ao 



'q^Ao 



a/3 



(20) 



(21) 



as shown in Appendix C. 

From ( p0| ) and ( pl| ) we see that Gaussian phase fluctuations do not affect transverse correlation functions. In 
particular, the superfluid stiffness is unrenormalized. However, longitudinal correlations are affected in general. 
While Ao does not satisfy the /—sum rule, A does, which implies restoration of gauge invariance. Further, from ( |2~l| ) 
we can see that x nas a pole for q — > unlike xo, which leads to a collective mode which we will discuss in the next 
Section. However, the static compressibility given by — %(q — * 0,u) n = 0) is unaffected at the Gaussian level. 

We clearly see that the gauge-invariant (q, ui) dependent correlation functions are different from the mean field 
correlations which appear as coefficients of the phase action. This is not surprising since the phase variable in fact 
contributes to the physical longitudinal correlation functions and modifies the mean field result. This is true even for 
charged systems as will be shown below. 



B. Charged systems 

In charged systems we have to take into account the long range Coulomb interaction with energy density 



G 



^-ffcoulomb — E VqPqP-q 

q 



(22) 



where p q = jj ff c]^ CT c k,CT is the electron density. In anisotropic layered systems is given by 16 3 



Vq = 



2-Ke 2 d c 
q\\£b 



sm 



cos 



h.(q\\d c ) - cos(qj_d c ) 



(23) 



where q\\,q± denote in-plane and c— axis components of q, d c denotes the mean interlayer spacing, and tb is the 
background dielectric constant. We assume qua <C 1 always, where the in-plane lattice spacing a = 1. For q\\d c , q±d c <C 
1 this reduces to the ordinary 3D result: Vq = Aire 2 / {tbq 2 ) ■ 



We take our Hamiltonian to be K + H D 



H, 



coulomb- 



Since the short range attraction of -ff pa ir is important for 



small center-of-mass momentum while the Coulomb effects are important for small momentum transfer, we believe 
the breakup of the actual interaction in this manner is physically sensible and does not lead to any "overcounting" . 
The Coulomb interaction can now be decoupled using a field Uq{r) as 



exp 



■ E VqWP-qW 

q>0 



V(Uq(T),Uq(T))eX P 



—UqUq + i ^2 ( W qP-q + W qPq) 



q>0 



(24) 



To integrate over independent modes, since U-q — we only sum over q > 0. The last term in (|24|) can be recast 
as ^ r U{v)p{v). Thus the density p couples to the scalar field U in the same way as it couples to d T 6 in (|l4|). 

On integrating out the fermions we arrive at an effective action for the phase 9 q and the scalar potential lA q , given 

by 



S[6,U] 



1 

8T 



E 

q,^Ti 



with 



M' 1 = 



[0*(q,»w») U*(q,iuj n )] M~ 



aj 2 nX o + A Q/3 q Q q /3 2iu nXo 



0(q, iw n ) 
U(q,iuj n ) 



(25) 



-2iu nXo * 4(- X o + Vq 1 ) 

where xo* = Xo(— q> —iw n ). Integrating out the field U leads to the action 

S charg e d [e] = ^ E (-^nX*^ + Ao^qV) 0(q, ^„)0(-q, -iw n ) 



(26) 



(27) 



q.oJn 



where the mean field charged system density correlator Xo = Xo/(l — KjXo)- This form of the phase action is 
independent of the order parameter symmetry and has been obtained earlier for s-wave SCs @ . We note that this 
form differs considerably from that assumed in Ref. |l5| ], where the physical (q, <x>)-dependent longitudinal dielectric 
function appears as a coefficient in the action. As emphasized earlier, physical longitudinal correlations cannot appear 
as coefficients of the phase action. 

The regime of validity for the above action is q < 7r/£ as for neutral systems (see discussion below (|l9|)). The 
frequency cutoff is given by \xo' PA (q)u 2 l \ < -Econd and thus depends on q. In particular, for q — > 0, the action remains 
valid even at frequencies larger than the gap and can be used to obtain the q — * plasma mode. 

The gauge-invariant correlations for the charged system are obtained from the neutral system results (|2(]) and ( |2~l| ) 
by replacing xo — * Xo > as discussed in Appendix C. 



VI. PLASMONS 



In the previous Section we have found phase actions of the form S[6] = ^ J2qu n 1 \Q( c L> UJ n )\ 2 for neutral sys- 
tems (see eq. (|l5|)) and for charged SCs (see eq. (|27|)). The dispersion of the collective phase mode is defined by 
TZeAlg 1 (q,u>) = 0. For neutral systems, we note that this condition is identical to demanding a pole in physical 
density-density correlation x given in (21). Phase and density fluctuations are thus coupled and share the pole of the 
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collective mode. This is true even for charged systems, where the plasma frequency, u> p , corresponds to the pole of 
the physical density correlator and is given by lim q ^o "^-ex~ 1 (q 1 w p ) = 0. 

The phase action is valid for all frequencies such that u> 2 < E con dVq and thus is valid even for large frequencies for 
q — * 0. If the plasmon is at finite frequency for q — ► 0, Landau singularities do not occur at finite temperatures. One 
can thus use this action to obtain the q — > plasma mode at zero and finite temperatures. Further, to have a sharp 
plasmon the damping must be relatively small: TmMg <C ui p . 

In this section, we first briefly consider neutral systems followed by a discussion of charged systems. For charged 
systems, we study the in-plane plasma mode and then consider c— axis plasmons for systems with a finite c— axis 
superfiuid stiffness. Our discussion of the c-axis plasmon is to a large extent independent of the details of any c— axis 
model. We also make an estimate of the c— axis plasma frequency for Bi2212 obtained within this phase action and 
compare it with experiment. 



A. Neutral Systems 

For neutral systems we have Mg 1 = |(Aq q^q^ — xo^n)- At T = 0, continuing to real frequency, the collective 
mode frequency obtained from above is given by uj(q) — cq where the sound velocity c = \J D a s j ' n. This reduces to 
the standard result c = v F /\/d in the weak coupling limit in the continuum, where d is the spatial dimension and v w 
is the Fermi velocity. At finite temperatures, there are Landau singularities in and xo a s seen from 
which prevents us from taking the q — > 0, lu — ► limit. 



B. Charged Systems 

In the action (|27|), we can set Xo FA ~ > — V^q m ^ ne nm it q — > 0. Analytically continuing to real frequency, and 
setting TZeMg 1 (c[,uj p ) = 0, we obtain 



lim 



0, (28) 



where uj p depends on the direction of propagation (in-plane or c— axis) in anisotropic systems. 

We first consider the in-plane plasmon. For q — > (q in-plane) and finite u> we can write = d a p {—iet, uja{uj)/e 2 ) 
where a = a' + ia" is the in-plane optical conductivity (including the effects of Gaussian phase fluctuations). The 
background dielectric constant enters in this definition of a since the conduction electrons are affected by the electric 
field which has been screened by the background. We thus see that 

to V ^J&^A = 47r ^ bK) = ^ (29 ) 
q— »o e 

<j{lo) thus governs the location and damping of the plasma mode through a self-consistent equation. This equation is 
identical to demanding a zero in the real part of the longitudinal dielectric constant (e), since e = e& + Amia /uj in a 
gauge invariant theory. 

At this stage it is instructive to compare (a) the physical plasma frequency u p , (b) the conductivity sum rule plasma 
frequency uj* defined by 

l rf ^'H = ^E^( k )W k )> s f' ( 3 °) 

and (c) the superfiuid plasma frequency cu ps defined by tu 2 s {T) = 4ne 2 D s (T)/et„ where D S (T) is the T-dependent 
superfiuid stiffness related to A(T), the penetration depth |0j . With this definition, the real part of a can be represented 
as o~'(uj, T) — (Lu 2 s (T)/4)S(uj) + a re g(uj, T) where o~ reg is the regular part. Finally, we have the Kramers-Kronig relation 
for a(uS): 

i /*°° o 
= ~V j Q *Jo>V )ljr -^ (31) 

We will now try to use these relations and the structure of <j'{ui) to obtain direct information about the behavior of 
the plasma mode, which is not directly seen in an optical conductivity measurement. 
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Conventional clean 3D s-wave SCs at T = have a very large superfluid stiffness which can be inferred from 
penetration depth measurements, and little spectral weight at higher energies. Ignoring interband transitions, we can 
then take cr'(uj) w (cj 2 . (01/4 Wo;) which leads to lo* — uj ps (0) from (|30|). It also implies (t"(lj) — w 2 s (0)/47ra; through 
the Kramers-Kronig relation (|3~t|). Using ( [?9|) we then get lu p = w ps (0), a large plasma frequency. In the presence 
of weak disorder and at finite T, uj ps decreases and the spectral weight in <j'{lo) redistributes, leading to finite o~(tu) 
over energy scales r t ~ , A which are the quasiparticle transport lifetime and SC gap respectively. Since lo p 3> itr j A 
to begin with, it is unaffected by this low energy redistribution. This is easy to see from d3l| ) above for <t"(cj) where 
we can set ui' w in the denominator for the region of interest, and this along with (|30| ) and (^9|) leads to ui p — to*, 
which is insensitive to the spectral weight redistribution. Further, the small o~'(w p ) implies a sharp plasmon in this 
case. Thus for conventional s-wave SCs, we finally arrive at to p ~ to* > tu ps (T). The last relation is satisfied as an 
equality only in a Galilean invariant system at T = 0. 

For the cuprate superconductors, the in-plane ct'(oj, T = 0) has the following features: (i) a condensate contribution 
(e 2 nD°(T — Q)/e b )6(uj) and (n) absorption by quasiparticles which has low frequency spectral weight (in a d- wave SC) 
with features around twice the maximum gap followed by other higher energy features. The condensate contribution 
along with the large low energy spectral weight coming from quasiparticles is expected to lead to the large plasma 
frequency, as in the case of conventional SCs above. Ignoring interband transitions in calculating lj*, we then arrive 
at uj p « to* > Lu ps . 

The normal state in-plane plasma frequency has been measured to be large (~ leV) in the cuprates Jl^| while the 
spectral weight rearrangement in a' in going from the normal state to the SC state is over smaller energy scales [fl8[ . 
The high energy normal state plasmon is thus expected to smoothly go over into a high energy SC state plasmon 
expected from our above discussion, similar to conventional SCs. 



C. Josephson Plasmons along c-axis 

To study the c— axis plasmon we assume a non-vanishing c— axis stiffness D ± . We set q along the c— axis in (|2£ 
leading to w PjC = Ait<j"{Lij PtC ). We know lo ps ^ c is small in the high T c systems as seen from the large c-axis penetration 
depths, and cr' c (uj) is measured to be very small over a very large energy range [[l9|. This is partly due to the form of the 



c-axis dispersion which is proportional to (cos k x — cos k y ) |20( ], so that nodal quasiparticles which lead to low energy 
spectral weight in-plane have a much smaller contribution along the c-axis. Thus at T = 0, the only "free carriers" 
come from the condensate, and to* ~ uj ps ^ c . On using ( |3l| ) and (|3(]), this leads to w PiC (T = 0) = LJ pStC (T = 0). With 
increasing temperature, the spectral weight in a" gets transferred to very high energies |21[ . The plasma frequency 
then continues to be given by w P]C (T) = uj ps ^ c (T), which decreases with increasing T and vanishes above T c as seen 
in experiment Jl^,^,^3|. Thus the c-axis plasma oscillations are seen only below T c , justifying the use of the term 
"Josephson plasmon" . 

A model which assumes disordered quasiparticle transport along the c-axis pi| and a model which only permits 
pair tunneling p5| both appear to lead to this behavior for uj pc (T). Presumably the main effect of disorder in the 
first model is to suppress single particle tunneling leading to pair tunneling as the dominant process below T c ; it then 
appears that the two models are similar in spirit. The disorder scale required to reproduce the experimental results 
in the first model appears to be large, of the scale of the one particle tunneling bandwidth. One reason for the large 
disorder scale could be that the c-axis dispersion in this calculation has been chosen to be independent of k x ,k y , while 
the actual (cosfc x — cosfc y ) 2 dependence would already suppress single-particle tunneling near the nodes in the clean 
case. This might then lead to a smaller disorder scale required to suppress this tunneling process completely. 

We finally proceed to compare the c-axis plasma frequency with the stiffness obtained from penetration depth 
experiments in Bi2212. We consider only the in-phase c-axis plasmon and ignore the "optical mode" corresponding to 
out-of-phase fluctuations arising from the bilayer structure, which is expected to be at a higher energy [ p6|p7| ]. This 
leads to 

9 47re 2 c 2 . , 

^=^=^1 (32) 

where c is the velocity of light and A c is the low temperature c— axis penetration depth. A c in Bi2212 has been 
measured p8| to be about 100 microns. Using this and setting ~ 10, we get uj p c w IK. This is in reasonable 
agreement with Lu PtC ps 8-10K extracted from experiment |29|-|3~tfl given experimental errors, and uncertainties in the 
estimate of e^. 







D. Plasmon dispersion 



In order to understand the plasmon dispersion and the variation of the plasma energy with direction of propagation, 
we consider a simplified model for the in-plane and out-of-plane conductivity. Since w p , a b ~ and is large for the 
in-plane plasmon, we set the in-plane conductivity <j'(u>) — (u>p /4)5(u>), consistent with the conductivity sum rule, 
which then reproduces the large in-plane plasma frequency. We emphasize that this simplified form for the in-plane 
conductivity is valid only at large frequency, for studying the high energy in-plane plasmon. It is not valid for 
considering low energy in-plane physical observables, such as the superfluid stiffness. Using this form for the real 
conductivity, cr" b (w) = w* 2 /47toj. The c-axis conductivity is given by a' c = {uj pa c /4)<5(u;) as discussed in the last 
Section, which leads to <j"(u>) = u) 2 s c /4iru). This simplified model is tailored to capture the correct plasma frequency 
for propagation along the in-plane and c-axis directions. For an arbitrary direction of propagation, we expect to 
obtain a reasonable interpolation. 

For general q, in the absence of detailed information on the cr(q, to), we assume that the conductivity is independent 
of q. This appears to be a reasonable assumption at high energy, and leads to the plasma frequency being given by 

w p(<0 = ^T (a': b (^ + a':^ p )ql) (33) 

The plasma frequency for |q| — > is a function of the angle of propagation, ip, measured with respect to the a6-plane 
and given by = (w* 2 cos 2 (ip) + w 2 sc sin 2 (?/>)) and varies smoothly from the low energy Josephson plasmon for 

c-axis propagation to the high energy plasmon in-plane. In order to examine the plasmon dispersion in a particular 
case, we consider the limit q ± = n/d c and study uj p as a function of (q,,cj). Normal state data for LSCO jl7| shows 
Wp,\\{q±. = 0,q y — ► 0) ~ leV; we assume a similar value for Bi2212 with ui Vt n(q± = 0,q y — > 0) w ui*. Using d c /a » 4 
and ujpc — LUps tC ~ l.OmeV p9|-|3l[|, we plot the plasma frequency in Figj^. This crosses over from a low energy 
Josephson plasmon for q» — > 0, corresponding to c-axis propagation, to a nearly two dimensional plasmon dispersing 
as ^/qj at larger q n . At small q^ , the dispersion is known to be acoustic for Lu pStC — fl32[ . It appears to be acoustic 
in Bi2212 at small q.. (see Fig. due the extremely small value of LOp S , c , but finally levels off leading to a finite 
Josephson plasmon gap for q^ — ► as shown in the inset in Fig.[j]. The 2D ^Jq^ dispersion is obtained mathematically 

in the limit of large d c /a and is given by lu p ,2D (q|| ) = ( W p/V2) y/q~» d c and we plot this in Fig.|l| for comparison. It must 
be emphasized that plasmon damping would be important in the real system and would have non-trivial dependence 
on the angle of propagation. The sharp plasmon we obtain in the above cases is an artifact of our simplified model 
for the conductivity. 



VII. QUANTUM XY MODEL 

The superfluid stiffness obtained above D°(T) = ^ [J2k m xx (1 ~ W-^k) — 2^ k v 2 (~df/dE k )] is unaffected by 
Gaussian phase fluctuations. Corrections to this result are unimportant in conventional superconductors which have a 
large coherence length and a large D®(T = 0). However, as we show below, effects beyond the Gaussian approximation 
could lead to large corrections in systems with a small coherence length and small £>!?. To study such effects we derive 
in this Section a quantum XY model and analyze quantum and thermal phase fluctuations in the following Section. 
For clarity of presentation, we outline the derivation for a (i-dimensional isotropic system; the generalization to the 
anisotropic case is straightforward. 

The quantum XY model describes the dynamics of the phase variables 0r(t) defined on a lattice with lattice spacing 
£o> the coherence length. The simplest action periodic under 0r(t) — * 0r(t) + 2ir is given by 

Sxy[6] = J2 A(Q)ujl\9(Q,tLu n )\ 2 +B dr £ [1 - cos(0r(t) - ^(r))] , (34) 

Q,o)„ <R.R'> 

where (R, R') are neighboring sites. It is important to emphasize that, given the cos(#r, — #r/) form, there are no 
constraints on the spatial gradient of the phases defined on the coarse-grained scale of £o- This is in contrast to the 
Gaussian action (^), derived on the scale of the microscopic lattice spacing a, which could only describe slow spatial 
fluctuations of the phase whose energy did not exceed the condensation energy. 

Our task now is to determine the coefficients A(Q) and B of this effective action. Unlike in some other cases 
p3| it is not possible to directly derive the quantum XY action from the underlying fermionic Hamiltonian, since 
the cumulant expansion we used to derive effective phase actions was controlled by the smallness of spatio-temporal 
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gradients in 0. We therefore proceed as follows: we compare the action ([m]) in the limit of slow spatial variations on 
the scale of the coherence length and match coefficients with those of the Gaussian action: 

S W = ^E"^!^'^)! 2 + I f 1/T drY,D° s a d - 2 %{r) - r+Q (r)] 2 (35) 

where V(q) is the generalized d-dimensional Coulomb interaction. For the (q, iu> n ) of interest, we have set Xq PA ~ 
— l/y(q) and Ao — D®(T). The double prime on the summation denotes momentum and frequency cutoffs 



|q| < <? c = tt/£o and u 2 n < (27m c T) 2 = D» [ - ) F q (36) 



which arise from demanding that the energy cost of the terms in (35) to be less than the condensation energy 
£-cond = |-Ds (^/Co) 2 - The q cutoff can also be viewed as a representation of the q-dependent stiffness being roughly 
constant (« D®) for q < q c and decreasing to zero for q > q c . 

In arriving at (|35|), we have Fourier transformed the gradient term from the (q, ito n ) to (r, t) variables. While the 
above r-local form of this term is true when all Matsubara frequencies are present, we now determine its regime of 
validity given the frequency cutoff in (^) . With the cutoff in ( |36| ) , the gradient term on Fourier transforming is given 
by 

| V D° s a d q 2 f 1/T dr dr'eiq, r)9(q, t')K(T(t - r')) (37) 

8 q JO 

In terms of the dimensionless quantity z = tT the kernel is given by K(z) — sin [(2n c + 1)ttz] / sin(7rz) where 
n c = n c (q) given by ([36|). The kernel K(z) is periodic in z, K(z + 1) = K(z), and it is sharply peaked around 
z = for large n c . The width of the peak can be estimated from the first zero of K(z) as z$ = l/(2n c + 1). Fra- 
nc £ 10, zq <C 1 which is true in the low temperature regime that we shall be interested in. We thus approximate 



K{z) as a delta function in "time" and work with a local-r action in (35). 

To determine the coupling B in (|3~4|), we consider static 9 configurations and apply a small external twist $ which 
will be distributed uniformly over the system. For the Gaussian model ([35]), with lattice spacing a — 1 and TV = L d 
sites, the phase twist per link is ($/L), while for the the XY model (p4[), on the "coarse grained" lattice with lattice 
spacing £ and (£/£o) d sites, there is a larger phase gradient <J>/(L/£o). Since the total energy cost for this phase 

twist is the same in the two cases, one obtains D°a d ~ 2 (®/ L ) L d /% = B {Ho/La) (La/£ ) /2, which leads to 

Similarly, for a r-dependent phase fluctuation at a frequency u> n , (L/£o) d phases contribute in the "coarse grained" 
XY model, as opposed to L d in the Gaussian case. We thus get A(Q)uj 2 [La/^o) d = (ui 2 a d /8Vq)L d on equating the 
first term in the two actions. This leads to A(Q) = (, d /8Vq. Finally, noting that the momentum Q = q£o/a since 
the distances in the "coarse grained" lattice are in units of £o, this can also be written as A(Q) — £q/8V"(Q), where 
V(Q) = V(Qo/fi)) = 7(«3). 

Having obtained A(Q) and B in the isotropic d-dimensional case, we generalize to anisotropic systems. We work 
with a layered 3D system with lattice spacing a = 1 in-plane and d c along the c-axis. The in-plane coherence length 
£o S> a and the c-axis coherence length £ x — d c . In this case, since £ x = d c , we proceed exactly as above but coarse 
grain only the in-plane variables. Denoting the Gaussian model stiffness by D® in-plane and D° ± along the c-axis, we 
arrive at the final action 



s w = w E wnM-Q, -<*o + ^r 1 t /T dT E (! - cos ^ t ) - ( r + a < T )i) ( 38 ) 

Q,w„ Q T,a=x,y 



D° x d c (toY f 1/T 



a 



f dTV(l-cos[0(r,T)-0(r + z,r)]) 
Jo 



where V(Q) = V(Q a/Co) and ^(Q) is the Coulomb interaction for layered systems given in (p3|). While all 
contribute in ( p9| ) above, the prime on the summation denotes the Matsubara cutoff consistent with 
couplings in this action depend crucially on £o and we shall examine the consequences of this below. 



Q| < 7T 

The 
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VIII. RENORMALIZATION OF THE STIFFNESS 



The quantum XY model action has both longitudinal fluctuations and transverse (vortex) excitations. Near a finite 
temperature phase transition, the dynamics is unimportant and we recover the classical XY model with the possibility 
of a phase transition in the 3D-XY universality class. In this Section, we examine low temperature properties and 
the effect of quantum dynamics. We ignore vortex-antivortex pair excitations since these have a core energy cost 
and would be exponentially suppressed at low temperatures. We deal with the longitudinal fluctuations within a self 
consistent harmonic approximation (SCHA). 

To examine the low temperature in-plane properties, we assume D a = in ( |39| ) since it is very small in highly 
anisotropic systems with a large X ± . The c-axis stiffness would become important if X ± < A. (£o/ a )i which implies 
D. < D ± (£o/a) 2 , leading to the c-axis contribution being important in (|39|). For Bi2212, detailed calculations, which 
we omit here, show that it does not affect our in-plane results. 

The SCHA [j34|]4| is carried out by replacing the above action by a trial harmonic theory with the renormalized 
stiffness chosen to minimize the free energy of the trial action. This leads to £) M = exp(— (S9 2 )/2) where 
SO = (#r iT — Or+a.r) and the expectation value is evaluated in the renormalized harmonic theory. Explicitly, we get 



2/ / d3 Q V ^ IVY) 

' ' ^)\^ n uie o d c /V Cl +D..d c e Q 



where e(Q) = 4 — 2cosQ x — 2cosQj / . 

We have analyzed the above equations to extract information about the importance of quantum and thermal phase 

fluctuations. Our numerical results can be simply summarized as follows: (S8 2 )(T = 0) ~ J (e 2 / £b^o) I (0)d c ) 

is a measure of quantum fluctuations, while thermal fluctuations become important above a crossover scale T x ~ 

(0)d c )(e 2 /tbio). These quantities are simply understood as the zero point spread and the energy level spacing 

of an oscillator respectively, in the renormalized harmonic theory. The scale T x is also the temperature at which 
n c(Q = 7r) ~ 1 (giving a broad kernel K(z)), with non-local effects in r becoming important. 

It is easy to see that phase fluctuation effects are negligible in the BCS limit of large £o- With e 2 /e^a ~ D., d c ~ E F 

and £ ~ v f/A, one obtains the standard result E cond ~ A 2 /E F per unit cell, and (59 2 )(T = 0) ~ ^ A/E F -C 1 and 
T x ~ ^E~K » T c . 

For the cuprates the short coherence length and small D„ act together to increase (SO 2 ), but they push T x in 
opposite directions. For optimal Bi2212 we use e 2 /eba w 0.3eV with tb ~ 10 and £o/a w 10. Bi2212 has a bilayer 
stacking structure with the planes within a bilayer being much closer than the distance between bilayers. Assuming 
the phase within the bilayer to be fully correlated, we set d c to be the mean inter- bilayer spacing and thus d c /a « 4. 
Using A,, (0) w 2100A, we then get the bilayer stiffness (0)d c ~ 75meV Q. These values lead to E con( i 6if/unitcell 
which is somewhat larger than estimates for optimally doped YBCO from specific heat measurements | fj5f ; we are 
unaware of similar data for Bi2212. We find that the crossover scale T x « 350K ^> T c . Since the bare stiffness 
D® actually decreases with temperature due to quasiparticle excitations, a better estimate of the thermal crossover 

scale may be obtained from T x ~ \J ^ x )( e2 / £ b£o)'-, this results in a crossover scale T x ~ T c . Thus, thermal 

fluctuations are clearly unimportant at low temperatures T <C T c . Quantum fluctuations are important since we find 
(S0 2 )(T = 0) ~ 1 at optimal doping. 

To study the temperature dependence of A„ (T) and the bilayer stiffness D ] .{T)d c ||, we set the bare stiffness 
d c D®(T) — d c D®(0) — 2a°T, where the linear decrease arises purely from nodal quasiparticle excitations within a 
single layer. This implies 1/A 2 (T) — 1/A 2 (0) — (47re 2 /h 2 c 2 d^j 2a°T. We plot the results of a numerical calculation 
of 1/A 2 (T) in Fig (||). Phase fluctuations are seen to lead to a large quantum renormalization of 1/A 2 (0) and to very 
little change in the slope of 1/A 2 (T) |3^]. The negligible renormalization of the slope of 1/A 2 (T) which we find, is true 
for a range of parameter values around our specific choice which has been constrained by experiments. It is however 
not the case more generally, and the slope could be renormalized by quantum fluctuations for a very different choice of 
parameter values. This effect of quantum fluctuations should be contrasted with the effect of classical thermal phase 
fluctuations which do not renormalize A y (0) or (0), but increase the slope of 1/A 2 (T) relative to its bare value. 

We note that in the absence of quasiparticles, the superfluid stiffness in this model would have an exponentially small 
temperature dependence, arising from phase fluctuations which are gapped. While it might appear that there could 
be low temperature crossovers resulting from the c-axis plasmon being at low energy (~ 10K for Bi2212), the phase 
space for these low lying fluctuations is extremely small to lead to a linear T behavior. Even in a purely 2D system, 
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which supports (gapless) low energy plasmons dispersing as y/gjj", the phase stiffness would decrease slowly, with a 
large power law (~ T 5 ). Quasiparticles are thus crucial in obtaining the observed linear temperature dependence. 
We find that we have to choose the bilayer stiffness d c D®(0) w 130meV corresponding to a bare A||,o(0) ~ 1600A 

and a slope a « 0.35meV/ K to obtain the renormalized values d c D^ ss 75meV (implying A,, ~ 2100A) and 
a ps 0.35meV/ K , in agreement with experiment. Thus quantum effects lead to a large (~ 40%) decrease of 1/A^(0), 
but no change in its linear T slope [|36f for our choice of parameter values. 

The bare slope a within a theory of non-interacting Bogolubov quasiparticles is given by (k B lrL2/ir)(hvFk F / Ad). 
Using the measured ARPES dispersion this leads to a ~ 0.8meV/K which is much larger than the 'bare' value 
we have used above to obtain agreement with penetration depth experiments. This points to the inadequacy of the 
non-interacting quasiparticle picture. The above discrepancy could be accounted for by considering quasiparticle 
interaction effects at the mean field level before considering the effect of phase fluctuations. These interaction effects 
become more important as one underdopes to approach the Mott insulator p7j. 



IX. CONCLUSIONS 



In this paper, we have focussed on the excitations of a short coherence length <i-wave superconductor. These are 
nodal fermions and the fluctuations of the amplitude and phase of the order parameter. Using an effective phase- 
only action we have discussed collective plasma modes and renormalization of superfluid stiffness by anharmonic 
longitudinal phase fluctuations. We summarize below some of our main conclusions. 

We have found that the important excitations are the low-lying fermionic states near the nodes and quantum phase 
fluctuations of the order parameter. Although the d-wave state supports in addition, two amplitude fields and a 
bond-phase field, these have been shown to have negligible spectral weight at low energy and are unimportant for 
the low temperature thermodynamics. They could possibly be probed in experiments such as Raman spectroscopy 
measurements designed to detect these fluctuations. 

Our discussion and derivation of the plasma modes emphasizes a unified way of looking at the in-plane and c— axis 
plasmons, in a manner which is relatively independent of detailed models of c— axis propagation. The very different 
nature of the two plasma modes, with a small c— axis plasma frequency governed by the c— axis stiffness and a large 
in-plane plasma frequency not directly related to the in-plane stiffness, can both be understood within our phase 
action. A microscopic derivation of c— axis conductivity sum-rules and T-dependent spectral weight transfers would 
depend on specific models |2^j2^], and we have not discussed these. 

Our derivation and treatment of the effective phase-only action emphasizes the crucial role played by the coherence 
length in imposing momentum and frequency cutoffs in the phase fluctuations. This allows us to interpolate from the 
BCS limit where phase fluctuations are unimportant to a regime of strong quantum fluctuations in the short coherence 
length limit. We find that quantum and thermal fluctuations cannot both be present at low temperatures; The short 
coherence length increases quantum fluctuations while pushing up the temperature scale at which one crosses over to 
thermal phase fluctuations. The strong longitudinal quantum fluctuations of the phase predicted by our calculation 
would also imply dynamical charge density fluctuations at low temperatures in the SC phase, which could possibly 
be probed in experiments. 

It has been pointed out that there is a discrepancy in the magnitude of the linear T slope of the measured penetration 
depth and the value calculated using ARPES data assuming free quasiparticles [0. However, phase fluctuations effects 
had not been taken into account before comparing data from the two measurements. We find that even including 
the relevant quantum phase fluctuations, a discrepancy is present which points to strong quasiparticle interactions 
even at optimal doping. A theory to account for these quasiparticle interactions is however lacking. One possibility 
is to invoke a phenomenological superfluid Fermi liquid theory description for the quasiparticles |58 
out however, that such a theory has a large number of free parameters and lacks predictive power although the 
experimental results may be easily rationalized. 

Finally, we have restricted our study in this paper to the low temperature properties of the SC without addressing 
the issue of what happens at higher temperatures within the SC state and in the normal state. This leads naturally 
to the problem of fermions interacting with a strongly fluctuating order parameter, which is at present an important 
open problem. 



Note added in proof: We have recently studied, in some detail, the effect of ohmic dissipation on phase fluctuations 
p9[ . Such dissipation, arising from a finite low frequency optical conductivity, is seen to reduce the magnitude of 
quantum fluctuations and reduce our estimate for the thermal crossover scale. Nevertheless, we find that the crossover 
scale is still large, so that our conclusion, about quasiparticles dominating the low temperature behavior of response 
functions, remains unchanged. 



13 



Acknowledgements: A. P. thanks D. Gaitonde and C. Panagopoulos for useful discussions. The work of M.R. 
was supported in part by the D.S.T., (Govt, of India) under the Swarnajayanti scheme. 



APPENDIX A 



In this Appendix, we present an approximate analysis of the density of states for the amplitude fields rj Std and the 
bond-phase field <f>. This analysis gives us insight into the nature of fermionic excitations which contribute to the low 
energy spectral weight for these fields and recovers the power law for the low energy DOS obtained in the numerics. 

A.l Amplitude fluctuations: 

The low energy density of states, Ni(u>) — jj ~ 2"mMj(q, uj)/tt i = s,d with the restriction |q x |, |q y | < tt/^q. 
The spectral weight XmJWi(q, u>) at T = arises from summing over low lying pair excitations at all momenta 
(k, k — q). This contributes to absorption at frequencies u> = E^ + -Ek-q- For q = 0, the spectral weight extends to 
uj = coming from low lying pair excitations from momenta k arbitrarily close to the node. At finite q the absorption 
sets in beyond a minimum threshold which corresponds to k at the node and k — q near the node (since q is small 



due to the momentum cutoff). The threshold is given by w m m(q) = y qfA 2 , + q^vp where q±, qi refer to components 

of q parallel and perpendicular to the Fermi surface respectively, at the node. 

We find numerically that the spectral weight for q ^ is nearly the same as for q = beyond the absorption 
threshold. We therefore approximate 

Ni(w) » ( — } [ A ° ^ImM,(q = 0, W )e( W - Wmra (q)) (40) 



where Q(x) is the unit step function which is 1 for x > and otherwise. For q = 0, the inverse propagator 

M7i - ^ - ^ry.Av , f* (4i) 



k 



where v?s(k) = cosfca, + cosk y and y>d(k) = cosfc^ — cosk y . Analytically continuing iu> n — > u> + i0 + and working at 
low frequency (w <C A<j) leads to 

- -TmM s4 = c 2 s A - w /2) ^ ( k ) (42) 

k 

where 2/c Sld = Adj? ^ k up d /E\^ — ip^ d ^/E^j. We evaluate the w-dependent momentum sum in ( |42| ) analytically 

by converting it to a Fermi surface integral and compute the constant c s-d numerically. Finally, doing the q sum to 
obtain the density of states leads to: 

N.(U>) _ =2 V 2 s {K.n) 



N qp (uo) S 167T A d v 



N d (uj) _ 2 1 



,4 



N qp {uj) d 256ttAX 

where ip s (k. Ftn ) refers to ip s evaluated at the gap node point on the Fermi surface and N qp (uj) = k F uj / {'Kv F lS. d ) is the 
quasiparticle DOS per spin which is linear in u>. We numerically estimate c s d ~ 10; The prefactors in the (E4) are 
then of order unity. 
A.2 The "bond-phase" field 0: 

For the 0— field, pair excitations similar to that for amplitude excitations lead to low energy spectral weight. This 
is easy to understand since both fields couple to the particle-particle channel with only different vertex factors. The 
behavior of ImM^,(q, w) is similar to that for amplitude fields with a vanishing threshold for q = and a finite 
threshold for q ^ 0. Following similar approximations, we set 

1 d 2 q 



N4u)= / -^ImM^(q = 0, W )e( W -w mm (q)) (44) 
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For T = 0, using particle-hole symmetry near the Fermi surface we find the inverse propagator 



M-( q =0, W „) = §^E 

k 



(45) 



Doing the integrals as before, we finally get 



n^w)/n„(u) = i- 4 ^(k F)B ) Jij- (46) 

where = ^ ^ k ^(k)/2J5k. The prefactor here is again of order unity. 

APPENDIX B 

In this Appendix, we briefly consider the linear time derivative term in the phase action, which we have dropped 
in the paper. For notational simplicity, we consider a neutral s-wave SC in 2D with lattice spacing a = 1. In carrying 
out the Hubbard Stratonovitch transformation as for the d-w&ve case, we introduce complex order parameter fields 
A r (r),A*(r) which are bosonic variables and satisfy the constraint A r (0) = A r (l/T). Writing A r = |A r |e r , this 
translates into |A r |(l/T) = |A r |(0), and 8 r (l/T) = 6 T (0) + 27rm r , and the partition function involves an additional 
sum over m r . Working in small 9 gradients and doing a cumulant expansion, one arrives at the following form of 
phase action for low momenta and frequencies: 



L dr ^ 

Jo 



i P e + \ne 2 + \d s (V6) 2 

o o 



(47) 



We now make the substitution 9(r,r) — 0(r,r) + 27rTm r r which implies 0(r, 1/T) = 0(r, 0). Substituting this in 
the action, we get 



S = 2ni p y"m r + ln 2 KTy"ml + ^y"(Vm r ) 2 + IttD s tS" [ ' dr r(V9) • (Vm r ) 

z — ' 2 z — ' bl £ — ' 2 In 

_ i- I- ' 



(48) 



+ J/ VT E(^ 2 +^( Ve ) 2 ) ( 49 ) 

where the derivatives denote discrete derivatives on the lattice. 

Now, at very low temperatures, D s /T ^> 1, we must set Vm r = while at high temperatures, kT ^> 1, we must set 
m r = 0. In either case, the field G decouples from the field m r and we get a Gaussian theory of phase fluctuations. 
The former condition (D s /T 3> 1) is equivalent to the condition that the spatial phase variation due to thermal effects 
is small; in particular, vortex configurations are unimportant. The latter condition (kT ^> 1) is just that the system 
starts behaving classically; since the extension along the imaginary time axis is 1/T, there is essentially no dynamics 
if 1/T — > and n is finite. 

In the presence of vortices, the core would described by a region where the magnitude of the order parameter |A r | 
decreases to zero. Since |A r | is a bosonic variable, the core would trace a closed loop in "time" 1/T. In this case, 
all the electrons which lie inside the loop undergo a phase change of 2ir each time the loop is traced while electrons 
outside the loop return to their original phase angle. This leads to a Berry phase factor iirpSiT) for the loop T with 
area S(T). The effect of this Berry phase factor on vortex dynamics was pointed out by Ao and Thouless [^T), except 
they obtained a coefficient of p s (the superfluid density) instead of p (the total electron density). Their result is 
special to a Galilean invariant systems at T = 0, where p s — p. Our result has also been derived earlier by Gaitonde 
and Ramakrishnan pu| . 

In the paper, we consider only slow spatial variations of the phase and work with just the periodic variable which 
we refer to as 9. We note that the linear time derivative term will not be important in the critical regime around 
the finite temperature superconductor to normal metal phase transition where dynamics is unimportant. It would 
however be important near quantum critical points at T = |42| . 
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APPENDIX C 



Physically relevant correlation functions should be gauge invariant. The functional integral method leads to a very 
simple and elegant way of demonstrating the role of phase fluctuations in restoring gauge invariance. (This is, of 
course, well known from early work of Anderson and others Jl3],[l4| ) . 

To this end we introduce external gauge potentials (A, Ao) in the Hamiltonian. This leads to the follow- 
ing modifications in the action (0). In the expression ( |l0| ) for Co we replace ^ by fj, + Aq(t, t) and 89/2 by 
89/2 — A rr /(r) where 89 = [0(r,r) — 0{r',r)]. Consider the gauge transformation Ao(r, t) — > Ao(r, r) + id T a(r, t) 
and A t] .i(t) — > A rr > (r) — (a(r', t) — a(r, r)). Invariance under this gauge transformation implies that the phase field 9 
must transform as 9(r, r) — > (9(r, r) + 2a(r, r) with | A| , <fi and the fermion fields unchanged. The correlation functions 
Xo and Ao"^ calculated at the mean field level, setting 9 = and <fi = are not gauge invariant. From the above 
discussion, it is clear that this problem can be solved by allowing for 9 fluctuations and integrating over these (rather 
than freezing 9 = 0). We now proceed to do this at the Gaussian level, which is entirely equivalent to the old RPA 
calculation. 

Physical correlation functions are obtained by integrating out the fermions and functional differentiation of the 
resulting Gaussian effective action with respect to the external sources Ao and A. We emphasize that these sources 
couple minimally to the original fermion operators, before the transformation C* — > c l e~ l9 / 2 in Section IV. We then 
find the density-density correlation allowing for Gaussian phase fluctuations to be 

qVAo^xo 
qVAo^-wfoo' 

This result obtains diagrammatically as follows. The fermions couple to the external source Ao and to the phase 
field 9 and we have to integrate out both 9 and the fermions. However 9 itself does not have a propagator unless the 
fermions are integrated out. To simplify the diagrammatic calculation we first introduce a fake term X9(q)9(—q) in 
the action which leads to the bare ^-propagator 1/A; we will take the limit A — > at the end. The diagrams in Fig.|| 
result for \- Summing the geometric series leads to 



X(q, iu n ) = z — af3 _ ^ . (50) 



X = Xo + 



Xo 2 ^l \ x _f Xo^l - A q 2 



A 



(51) 



(52) 



Taking the limit A — > then leads to the result in (|50|). 

We note that the physical static compressibility is given by x(q - * 0,cj„ = 0) = Xo(q - * 0, w n = 0), i.e., the 
field result is unaffected by phase fluctuations at the RPA level. 

Similarly, denoting the physical (q, a;) dependent stiffness by A, we get 

p a0 AcTAo-Vq" 
+ K Xo - Ao^q/'q"] 

We see that the transverse phase stiffness, (for instance along the a;— direction), given by A xx (cu n = 0, qj_ — ► 0, q x = 0) 
is unaffected by the Gaussian phase fluctuations. However, the longitudinal part of the current correlation function 
is affected, and 

~{jx3x)(u n = 0, q± = 0, q x 0) = ^=(j x j x )a(uJ n = 0, q ± = 0, q x 0) + A xx (o;„ = 0,q x = 0, q x 0) (53) 

k 

now satisfies the /-sum rule (which was violated at mean field level). Thus gauge invariance is restored. 

The derivation remains unchanged in charged systems, with the only difference being that we have to make the 
replacement xo ~ * Xq' PA = Xo/(l ~ KjXo) i n a H the equations (in this Appendix). This can be easily understood by 
comparing (|l5|) and (|27|) in the text. It is not hard to show that the longitudinal conductivity (o~l) defined through 
the longitudinal dielectric function by 

f ' = 1 IT/ = 1 + ( 55 ) 

in a gauge invariant theory, and the transverse conductivity defined by ctt(w) = iAe 2 / (tb^) (with A being the 
transverse part of K a/3 ) are equal in the limit q — > 0. In the text, we omit subscripts and refer to both conductivities 
by a since we work at q — > 0. It is easy to see that o~t{u) and hence o~(uS) is unaffected by phase fluctuations within 
RPA. However, it is only in a gauge invariant theory, such as the RPA, that (|5^) holds since it obtains from using the 
current conservation equation. 
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12 3 

FIG. 1. The plasmon dispersion in a layered system (Bi2212), as a function of q (| with q ± — n/d c (see text for details and 
parameters used). The ^/q^ dispersion in the 2D limit is plotted for comparison. The dispersion appears nearly acoustic for 
small q due to the very small Josephson plasma frequency (~ 10K) but rapidly crosses over to high energies (~ eV) with 
increasing q.. . Inset shows the behavior of the dispersion (in meV) for q.. — » 0. 
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FIG. 2. The bare and renormalized 1/A^ plotted as a function of temperature near optimal doping for Bi2212. We have 

chosen bare values such that the renormalized results, A y (0) w 2100A and dAy /dT « IQ.0A/K, are in approximate agreement 
with experiment. 




FIG. 3. Low order diagrams for \: The wavy lines indicate Ao, the external scalar potential. The dashed lines contribute 
1/A, the "fake" ^-propagator. The heavy dots refer to vertex factors of tu n arising from the vertex ipd T 6. Finally, the bubble 
contributions arising out of integrating out fermions are explicitly indicated. 
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